T1ρ, T2 and T2* mapping of lumbar intervertebral disc degeneration: a comparison study

Background Early and accurate assessment of lumbar intervertebral disc degeneration (IVDD) is very important to therapeutic strategy. This study aims to correlate and compare the performances of T1ρ, T2 and T2* mapping for Pfirrmann grades and morphologic changes in the IVDD. Methods This prospective study included 39 subjects with 195 lumbar discs. T1ρ, T2 and T2* mapping were performed, and T1ρ, T2 and T2* values of nucleus pulposus (NP), and anterior and posterior annulus fibrosus were measured. IVDD was assessed with Pfirrmann grading and morphologic changes (normal, bulging, herniation and annular fissure). The performances of T1ρ, T2 and T2* relaxation times were compared for detecting early (Pfirrmann grade II-III) and advanced degeneration (Pfirrmann grade IV–V), as well as for morphologic changes. Results T2 relaxation times was strongly corelated with T1ρ and T2* relaxation times. Areas under the curves (AUCs) of T1ρ, T2 and T2* relaxation times of NP were 0.70, 0.87 and 0.80 for early degeneration, and 0.91, 0.95 and 0.82 for advanced degeneration, respectively. AUCs of T1ρ, T2 and T2* relaxation times of NP were 0.78, 0.83 and 0.64 for bulging discs, 0.87, 0.89 and 0.69 for herniated discs, and 0.79, 0.82 and 0.69 for annular tearing, respectively. The AUC of T2 relaxation time was significantly higher than those of T1ρ relaxation times (both P < 0.01) for early IVDD, and the AUCs of T1ρ and T2 relaxation times for assessing advanced degeneration and morphologic changes were similar (P > 0.05) but significantly higher than that of T2*relaxation time (P < 0.01). Conclusions T2 mapping performed better than T1ρ mapping for the detection of early IVDD. T1ρ and T2 mapping performed similarly but better than T2* mapping for advanced degeneration and morphologic changes of IVDD. Supplementary Information The online version contains supplementary material available at 10.1186/s12891-022-06040-y.


Introduction
Low back pain affects up to 40%-80% of the population at some point during their lifetime, with considerable negative impacts on quality of life and social economy [1,2]. Intervertebral disk degeneration (IVDD) is considered to be one of the primary causes of low back pain [1,3]. Early phase of IVDD presents with biochemical changes, including decreases of overall proteoglycan and water content in the nucleus pulposus (NP). The later stages of IVDD are manifested in morphological changes, including disc bulging or herniation, annulus fibrosus (AF) tears, and vertebral osteophytes [4]. Unfortunately, current surgical treatments such as spinal fusion and replacement of disc are limited to pain relief treatment for severe IVDD. Surprisingly, recent studies showed regenerative strategies such as biologic and cell therapies may be applied for early-stage IVDD to limit disease progression [5,6]. In this context, noninvasive imaging Yang et al. BMC Musculoskeletal Disorders (2022) 23:1135 modalities are required for the detection of early changes and the severity of IVDD.
MRI is an important tool for evaluating IVDD. The Pfirrmann semi-quantitative classification is the most widely used for visual grading of IVDD [7], which is focused on the signal intensity and homogeneity of NP, distinction of the NP and AF, and height of the disc on sagittal T2-weighted images [8]. However, this scoring system is subjective, and it is unable to evaluate early biochemical changes of IVDD [9]. Besides the morphologic changes of disc bulging and herniation, high-intensity zone (HIZ) in the posterior annulus fibrosus (PAF) on T2-weighted images has been proposed to be an important sign for annular disruption and low back pain [10].
Since T1ρ, T2 and T2* mapping have different mechanisms, thus they may have different sensitivity to compositional changes throughout IVDD [9], we hypothesized that these MRI parameters are fundamentally different when assessing disc health. However, no studies have been performed to evaluate the relative performance of T1ρ, T2 and T2* relaxation times for both lumbar disc degenerative grades and morphologic changes. The purpose of this study was to assess the correlation and compare the performance of T1ρ, T2 and T2* mapping for Pfirrmann grades and morphologic changes in the lumbar IVDD, thus to determine the most appropriate method for lumbar IVDD assessment.

Study population
The study was approved by the institute's research ethics committee, and written informed consent was obtained from each participant prior to enrollment. Subjects with body mass index > 25, spinal fractures, tumors, infections, metabolic disease, previous lumbar disc surgery or interventional treatment, and contraindications for MR imaging were excluded. This study included 39 subjects (male 27, female 12; median age 44.0 years, age range 22-80 years) with single or recurrent episodes of low back pain and excluding other spine diseases except IVDD.

MRI examinations
All subjects underwent MRI in the morning to reduce the potential impact of diurnal changes in the discs. Lumbar MRI including sagittal T2-weighted imaging, T1ρ mapping, T2 mapping, T2* mapping and axial T2-weighted imaging were performed at 3.0 T scanner (Achieva, Philips Healthcare, Best, the Netherlands) with a dedicated 15-channel SENSE spine coil. The imaging parameters are shown in Table 1. T1ρ-weighted images were obtained with a rotary echo spin-lock pulse embedded in a three-dimension balanced fast field echo sequence. Spin-lock frequency was set as 500 Hz, spin lock durations were 0, 10, 20, 30, 40 ms.

Image analysis
A total of 195 discs covering L1-L2 to L5-S1 were analyzed. The morphological analyses were performed by the consensus of 2 experienced radiologist without knowledge of the quantitative MRI measurements. The assessments included Pfirrmann classifications on sagittal T2-weighted images [8], as well as the presence of disc bulging, herniation (protrusion or extrusion) according to the Lumbar Disc Nomenclature 2.0 [27], and annular fissure with HIZ of PAF on sagittal and axial T2-weighted images according to Aprill and Bogduk [9]. For Pfirrmann grades, IVDD with Pfirrmann grade I was normal disc, IVDD with Pfirrmann grade II-III were labeled as early degeneration, and IVDD with Pfirrmann grades IV and V were labeled as advanced degeneration [23].
The raw data of T1ρ mapping were fitted on a pixel-bypixel basis to the exponentially decaying T1ρ function using IDL 6.3 (ITT Visual Information Solutions, Boulder, CO) to generate a T1ρ map, and the T2 and T2* maps were generated and analyzed by using Image J software (National Institutes of Health, Bethesda, MD). Regions of interest (ROIs) were manually placed on the NP, anterior annulus fibrosus (AAF) and PAF of the T1ρ, T2, and T2* maps (Fig. 1). The anterior and posterior 20% of the disk diameter were labeled as AF with the remaining central 60% representing the NP [23,28]. Values of AAF and PAF were averaged as the value for the AF. All the ROI measurements were performed twice by an observer with an interval of one month.

Statistical analysis
Numerical data were expressed as median and range as the data was not normally distributed according to the Shapiro-Wilk test. For the reliability of two measurements of T1ρ, T2 and T2* relaxation times, the intraclass correlation coefficients (ICCs) were calculated, and ICC values that were ≥ 0.75 were considered excellent agreement [29]. Wilcoxon test was used to compare the quantitative MRI parameters between NP and AF, as well as between discs with or without HIZ in the PAF. Pearson correlations were used to evaluate the relationship between every two quantitative MRI parameters, and between quantitative MRI parameters and subject age. Spearman rank correlations were applied to assess the association between quantitative MRI parameters and Pfirrmann grades, or morphology, as well as between  subject age and Pfirrmann grades. The differences of quantitative MRI parameters in various Pfirrmann grades and morphology were evaluated with Kruskal-Wallis test and post-hoc tests. Furthermore, receiver operating characteristic (ROC) curves of each quantitative MRI parameter were plotted for detecting early IVDD (Pfirrmann grade II-III) and advanced IVDD (Pfirrmann grade IV-V), as well as for disc bulging or herniation, and annular tearing, and the cutoff values were defined with their sensitivities and specificities. The areas under the curves (AUCs) with 95% confidence intervals (CIs), sensitivities, and specificities were used to assess the diagnostic performances of quantitative MRI parameters. Delong method was used to compare the AUCs [30]. All statistical analyses were conducted using SPSS software (v23.0, IBM, Chicago, USA) and Medcalc (v18.2.1, Mariakerke, Belgium). A P value < 0.05 was regarded to be significantly different, and a Bonferroni-adjusted test was used for multiple comparisons.

Measurements of T1ρ, T2 and T2* relaxation time
The reproducibility between the two measurements for T1ρ, T2 and T2* relaxation times at different ROIs was excellent (ICCs = 0.88 ~ 0.93 T1ρ, T2 and T2* relaxation times of NP and AF were summarized in Table 2 and Table E4, and illustrated in   Fig. 3. The T1ρ, T2 and T2* relaxation times of NP significantly decreased with the increase of Pfirrmann grades (rho = -0.73, P < 0.001; rho = -0.88, P < 0.001, rho = -0.61, P < 0.001; respectively). Regarding T1ρ relaxation times of NP, Kruskal-Wallis test and post-hoc tests showed significant differences between all pair-wise comparisons of Pfirrmann grades (all P < 0.005) except for between grades I and II and between grades IV and V (both P > 0.005). For T2 relaxation times in NP, Kruskal-Wallis test and post-hoc tests showed significant differences between all of the Pfrrmann grades (all P < 0.005) except between grades IV and V (P > 0.005). Kruskal-Wallis test and post-hoc tests of the T2* relaxation time of NP showed significant differences between all pair-wise comparisons of Pfirrmann grades (all P < 0.005) except for between III and V and between grades IV and V (both P > 0.05).
The ROC curves of T1ρ, T2 and T2* relaxation times for early and advanced IVDD are presented in Fig. 4 and the corresponding diagnostic performance are provided in Table 3 and Table E5. The AUCs of T1ρ, T2 and T2* relaxation times of NP were 0.70 (95%CI: 0.59-0.80), 0.87 (95%CI: 0.80-0.93) and 0.80 (95%CI: 0.70-0.89) for detecting early IVDD, and 0.91 (95%CI: 0.87-0.95), 0.95 (95%CI: 0.92-0.98) and 0.82 (95%CI: 0.75-0.88) for identifying advanced IVDD, respectively. For assessing early IVDD, the AUC of NP in T2 value was significantly higher than that in T1ρ value (P = 0.001), but the AUCs of NP in T2 and T1ρ value were not significantly different from that in T2* values (adjusted P > 0.05). For assessing  IV-V, B). Numbers are areas under the curves with 95% confidence intervals in parentheses advanced IVDD, the AUC of NP in T2 value was significantly higher than that in T2* values (P < 0.001), but the AUCs of T2 and T1ρ values did not showed significant difference (adjusted P > 0.05). The AUC of NP in T1ρ value was significantly higher than that in T2* values (P = 0.004).

T1ρ, T2* and T2 relaxation times for disc bulging, herniation and annular tears
Of all 195 discs, 50 were labeled as bulging discs, 21 were herniated discs; 45 showed an HIZ in the PAF, reflecting annular tearing. T1ρ, T2 and T2* relaxation times of NP and AF for IVDD morphologic changes were showed in Table 4 and Table E6, and presented in Fig. 5. Kruskal-Wallis test and post-hoc tests showed that the T1ρ, T2 and T2* relaxation times of NP in normal discs were significantly higher than those in bulging or herniated discs and IVDD with HIZs in the PAF (P < 0.01), while there were no significant differences between bulging discs and herniated discs in all the MRI parameters (P > 0.05).
The ROC curves of T1ρ, T2 and T2* relaxation times for morphologic changes of IVDD are presented in Fig. 6, and the corresponding diagnostic performance are provided in Table 5 and Table E7. Calculating a ROC for discriminating bulging discs from normal discs, the AUCs of T1ρ, T2 and T2* relaxation times of NP were 0.78 (95%CI: 0.71-0.85), 0.83 (95%CI: 0.76-0.90) and 0.64 (95%CI: 0.55-0.72), respectively. Moreover, the AUCs of T2 relaxation times was significantly higher than that of T1ρ relaxation time (P < 0.01), and the AUCs of T1ρ and T2 relaxation times were significantly higher than that of T2* relaxation time (P < 0.01).
Calculating a ROC for identifying HIZ of PAF from discs without tear, the AUC of T1ρ, T2 and T2* Table 4 T1ρ, T2 and T2* relaxation times in nucleus pulpous for lumbar disc bulging, herniation and annular tears Values are expressed as medians and ranges in parentheses a Spearman rank correlations were applied to assess the association between quantitative MRI parameters and morphologic changes b Mann-Whitney U test was used to compare discs with and without annular tears

Discussion
Our study demonstrated significantly moderate to strong negative correlations between the three MR quantitative parameters (T1ρ, T2 and T2* relaxation times) of NP and Pfirrmann grades, and T2 relaxation times was strongly corelated with T1ρ and T2* relaxation times.
For identifying early IVDD (grade II-III), T2 mapping yielded greatest AUCs with high sensitivities, followed by T2* mapping and T1ρ mapping, indicating T2 mapping may be more accurate for the early disc degenerative changes. Additionally, T1ρ, T2 and T2* values of NP in normal discs differed significantly from those in bulging disc, herniated disc and annular tears. The AUCs of T1ρ and T2 mapping of NP were significantly higher than that of T2* mapping in detecting advanced IVDD (grade IV-V), bulging discs, herniated discs and annular tearing, implying T1ρ mapping and T2 mapping are more accurate than T2* mapping for late changes during IVDD. Our present finding is consistent with previous studies that T1ρ, T2 and T2* relaxation times of NP decreased linearly with increasing Pfirrmann grades [11][12][13][21][22][23]. Histologically, T1ρ is strong affinity with proteoglycan content in the disc matrix [13], T2 relaxation time has strong positive correlations with water and glycosaminoglycans content [19], and T2* relaxation time is highly correlated with glycosaminoglycan content [20]. During the IVDD process, the proteoglycan quantity diminishes and water content decreases [26], along with significant drop of T1ρ, T2 and T2* relaxation times of NP accordingly. T1ρ relaxation time of NP at Pfirrmann grade I was significantly higher than that at Pfirrmann grade III, IV and V, but there was no significant difference between grade I and II. However, there was significant differences in T2 and T2* relaxation times of NP between all pair-wise comparisons of Pfirrmann grades except for between grade IV and V for T2 and T2* mapping and between III and V for T2* mapping. The AUC of NP in T2 value was significantly higher than that in T1ρ value for assessing Pfirrmann grade II-III IVDD, with the cut-off values being similar to the Nagy et al. study [31]. This representation indicates T2-mapping may outperform the T1ρ mapping for detecting early stage IVDD. Yoon et al. reported T2 relaxation rates had stronger correlation with Pfirrmann grades than T1ρ value [9]. Menezes-Reis et al. [32] proposed that T2 mapping may be more appropriate than T1ρ mapping for the detection of early disc aging changes. Additionally, an ex vivo study showed T2 mapping was more sensitive to early degenerative changes than T1ρ mapping [19].
With regard to advanced degeneration and morphologic changes of IVDD (bulging discs, herniated discs and HIZs in the PAFs), the AUCs of NP in T2 and T1ρ mapping were similar (AUCs: 0.78-0.95) but significantly higher than that in T2* mapping (AUCs: 0.64-0.82). With ongoing IVDD, hydrophilic glycosaminoglycans within the NP is reduced or absent, result in structural defects such as annular tears and disc herniation [4]. Therefore, T2 and T1ρ mapping appeared to be sensitive to not only proteoglycan and water content but also collagen integrity. In line with this statement, a cadaveric study showed that the biomechanical properties of lumbar intervertebral discs are correlated with the collagen structure integrity [33].
Our results showed that T2 relaxation times was strongly correlated with T1ρ and T2* relaxation times. This may imply T2, T1ρ and T2* relaxation times are correlated to each other fundamentally, as all of them are associated with proteoglycan and water content [13,19,20]. A previous report also showed a significant linear association between T2 and T1ρ relaxation rates [9]. The results from our study suggest T2 mapping could be implemented as a synergistic modality with conventional MRI for assessing the severity of IVDD, which may help treatment option and prognosis evaluation of IVDD.
There are some limitations in present study. First, the number of patients included in this study was relatively small. However, power analysis depicted a minimum required number of 35 subjects with an alpha value of 0.05 and a specificity of 90%. Second, we only used Pfirrmann grade to assess IVDD, this grading system is subjective and cannot fully reflect the composition change during IVDD [33]. Meanwhile, histologic or biochemical validation of IVDD was lacking because we did not obtain specimens from the subjects. Further studies with humans and cadavers to validate the composition associated with different IVDD are warranted. Third, the echo times used in T2 mapping are relatively low, which may result in noise floor in T2 map, however the cut-off values of T2 mapping for early and advanced IVDD were similar to the Nagy et al. study [32]. Fourth, we did not correlate the T1ρ, T2 and T2* relaxation times with clinical symptom. Further study with a large cohort including clinical assessment is needed to clarify the relationship between MR quantitative parameters and clinical symptom.

Conclusions
In conclusion, our resulted showed T2 mapping performed better than T1ρ mapping for the detection of early IVDD, and T1ρ and T2 mapping performed similarly but better than T2* mapping for advanced degeneration and morphologic changes of IVDD. T2 mapping may be of great utility for detecting the early and later changes of IVDD, as well as for monitoring the treatment response to emerging regeneration therapies.
Additional file 1: Table E1. Correlation between subject age and T1ρ, T2* and T2 relaxation times of NP for each spinal level. Table E2. Correlation between age and T1ρ, T2* and T2 relaxation times of AAF for each spinal level. Table E3. Correlation between age and T1ρ, T2* and T2 relaxation times of PAF for each spinal level. Table E4. T1ρ, T2 and T2* relaxation times in AAF, PAF and AF for different Pfirrmann grades. Table E5. T1ρ, T2* and T2 relaxation times in AAF, PAF and AF for different Pfirrmann grades. Table E6. T1ρ, T2 and T2* relaxation times in AAF, PAF and AF for lumbar disc bulging, herniation and annular tears. Table E7. Diagnostic performance of T1ρ, T2 and T2* relaxation times in annulus fibrosus for lumbar disc bulging, herniation and annular tears.